Relativity accommodates superluminal mean velocities 
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Contrary to a widespread belief, measures of velocity can yield a value larger than c, the instan- 
taneous light speed in vacuum, without contradicting Einstein's relativity. Nevertheless, the effect 
turns out to be too small to explain the recently claimed superluminal velocity by the OPERA 
collaboration. Several other general relativistic effects acting on the OPERA neutrinos are also 
analyzed. All of them are unable to explain the OPERA result. 
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The recent observation of superluminal neutrinos by 
the OPERA team [1] (see also the MINOS collabora- 
tion [2]) has prompted an intense discussion. Particle 
physics arguments [3,4] or astrophysics data [5,6] seem 
to indicate that neutrinos could hardly attain such ve- 
locities. Moreover, and generally speaking, it is accepted 
that within the present physics knowledge, such veloci- 
ties cannot be contemplated if the validity of the basic 
laws of relativity is to be maintained. In particular the 
causal structure of spacetime would be badly violated if 
superluminal particles existed. 

However it should not be forgotten that such a causal 
structure holds good in special relativity, namely, in an 
ideal world where gravity is absent. When instead grav- 
ity is taken into account, a special relativistic spacetime 
arises only locally, in a close neighbourhood of every 
spacetime event P. Strictly speaking, with infinite math- 
ematial precision, this is true exactly at P (although, 
considering the inevitable inaccuracy of measurement de- 
vices, the region of validity of such an approximation 
turns out to be more or less extended, depending on the 
magnitude of the spacetime curvature at P). Thus, it 
should not come as a surprise that measurements carried 
out in a non-local way may apparently violate the main 
tenets of special relativity. 

We are going to show that, although the instantaneous 
velocity of a particle cannot be larger than c, relativity 
consents that the mean velocity of the same particle may 
well exceed c. The interesting and appealing aspect of the 
above statement is that what the OPERA collaboration 
has measured is precisely the neutrino mean velocity. 

To illustrate the above effect it will be assumed that 
the spacetime is satisfactorily described by the metric 
ds^ = gaodt^ — Qijdx^dx^ where {x^,x'^,x^) are gen- 
eral spatial coordinates and t is the coordinate time. 
Throughout the paper mass and time units will be cho- 
sen such that c and the Newton constant G are 1. Also, 
the Einstein summation convention for pairs of repeated 
indices will be used with Latin indices indicating spa- 
tial components and Greek indices referring indistinctly 
to both spatial and time components. It will be also 
supposed that all metric components are time indepen- 
dent, dgfj^n/dt — 0, and that the time-space ones vanish, 
g^i = 0. This last condition is not essential (we include 



it to avoid lengthy mathematical expressions) and can 
be omitted without changing the conclusions of the pa- 
per [7]. 

Consider a massless (or nearly massless) particle 
travelling along a spatial trajectory parametrized by 
a variable A and described by the three functions 
(x^(A),x^(A),x'^(A)). It travels from the spatial point Pi 
(where the parameter takes the value Ai) to the spatial 
point P2 (value A2). After having measured the physical 
distance A£ between the two points (the metric is time 
independent), an observer at P2 registers the time Ar 
used by the particle to complete the trip. By physical 
distance we understand what one obtains for instance by 
lying rods in succession between Pi and P2 (a more re- 
alistic procedure to extract A£ will be described in the 
discussion of the OPERA experiment). The time reading 
requires a previous synchronization among the clocks at 
Pi and P2. Then, the observer at P2 defines the mean 
velocity of the particle as 

V ^ M 
c~ At ' 

Let us see what general relativity predicts for this quan- 
tity. 

The physical distance A£ is given by 

A£ = y dX^Jgijx'i^ , (2) 



(1) 



where = dx^/dX. Viewed as a quadratic form, gij is 
positive definite. On the other hand, one of the several 
equations that govern the motion of the particle stems 
from the nullification of its proper time, ds^ = = 
goodt^ — g.ijdx'^dx^ , whence 



At = 



dX 



Vdijx'x^ 



(3) 



Expression (3) gives the coordinate time interval em- 
ployed by the particle during the trip from Pi to P2. 
To recover the time reading At on the observer's clock 
at P2 one has to multiply At = y/goo{X2)At. Therefore 
the mean velocity recorded by the observer is 
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The dependence of in (2)-(4) on A comes through 
their dependence on the spatial coordinates. 

The relevant fact that we wish to emphasize is that the 
mathematical framework of general relativity does not 
contain any mechanism whatsoever able to oblige (4) to 
stay equal or less than 1. Rather, as remarked above, (4) 
tends to 1 only when the velocity measurement becomes 
local, to wit, when po ints Pi and P2 tend to coincide. In- 
deed, since ^gijx'^xJ is positive, the mean value theorem 
can be applied and it states that a value Aj, comprised be- 
tween Ai and A2 , exists such that J^^ dX \/ g%jX^xi j gm = 

(l/-\/5oo(Aj)) dX gijX^xK Inserting this expression 
in the denominator of (4) we obtain 



c a/5oo(A2) 



(5) 



which tends to unity as Pi P2 because in this limit 
A J — A2 also. 

Let us consider an instance of the effect just ex- 
posed by studying the radial motion of a particle in a 
Schwarzschild metric (it characterizes the vacuum out- 
side a spherically symmetric mass distribution with to- 
tal mass m) in standard coordinates, ds"^ = goodt'^ — 
Qrrdr^ — r'^dfl^ where dft is the solid angle and goo = 
l/ffrr = 1 — 2m/r. The massless or almost massless par- 
ticle will travel radially from ri to r2 (we do not specify 
yet whether ri < r2 or r2 < ri). Because goo increases 
with r, following (5), we deduce u < c if r2 > ri and 
u > c if r2 < ri. The detachment from the light speed 
can be better appreciated by plotting v/ c against ri and 
r2- Specifically, for the Schwarzschild metric we have 
(calling = r^/(2m), A =1,2) 
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The absolute values in (6) make these expressions valid 

for both cases, r2 > ri and r2 < ri. This is an important 
remark because the fact that the time intervals for the 
particle to go from ri to r2 or for coming from r2 to ri 
are the same allows to reliably synchronize clocks with 
the exchange of light signals. 

In Figs. 1 and 2 the mean velocity w/c is plotted against 
1/pi for various values of p2- Note (i) that we have delib- 
erately excluded the region 1/pi > 1 because otherwise 
the particle would enter the horizon of the corresponding 
black hole and (ii) that all plots touch the line v/c = 1 
when pi coincides with p2 , confirming the validity of spe- 
cial relativity at short distances. 

The case ri > r2 (the particle approaching the spheri- 
cal distribution of mass) is shown in Fig 1. As viewed by 



the observer in r2 , w turns out to be always larger than c. 
This effect is particularly pronounced when the observer 
is close to the Schwarzschild radius = 2m (lower values 
of P2) and for large separations ri — r2. For p2 > 100 the 
resulting v/c is so close to 1 that the related plots cannot 
be seen within the scale of the vertical axis. 

In Fig. 2 the case ri < r2 (the particle going away from 
the spherical massive object) is treated. Now all mean 
velocities come smaller than c, which looks quite surpris- 
ing for massless particles. Again v tends to c whenever 
the radial coordinate ri of the position of the particle's 
departure is close to the radial coordinate r2 of the posi- 
tion of arrival. Note the tendency of v to become null as 
ri approaches for every r2- 
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FIG. 1. Mean velocity v/c oi the particle travelling from 
ri = 2mpi to r2 = 2mp2 in a Schwarzschild metric for several 
values of p2 as a function of 1/pi and always with n > r2. 
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FIG. 2. Mean velocity v/c of the particle travelling from 

ri = 2mpi to r2 = 2mp2 in a Schwarzschild metric for several 
values of p2 as a function of 1/pi and always with r2 > ri. 

As shown in the above plots, the supcrluminality or 
subluminality of the particle's mean velocity is dramat- 
ically enhanced in regions of large curvature (near the 
Schwarzschild radius r^). This stresses the general rela- 
tivistic character of the the described effect. 

A terrestrial experiment prepared to detect values for 
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the radial velocity v different from c cannot approach 
any Schwarzschild radius. If such an experimental set 
up was constrained to use ri, r2 > r©, the Earth's radius, 
then A£ and At in (6) could be approximated with excel- 
lent accuracy to linear order in m© / (m© is the Earth 
mass) and (1) would barely differ from 1, 



1 



m© 



TO© 



n - ; 



■ log — > 1 



(7) 



where the observer has been put on the Earth's surface, 
r2 = r©, and we have taken n > r©. For n not much 
larger than r©, the percentage of excess of velocity is a 
meagre {v — c)/c ^ 10^^°. 

To exemplify the above findings, we apply them in two 
cases: first, to the already mentioned OPERA experi- 
ment and, second, to the determination of the distance 
from the Earth to the Moon by exchange of light pulses. 

In the OPERA experiment a beam of muon neutri- 
nos was produced at CERN SPS and sent to the Gran 
Sasso underground laboratories (LNGS) in Italy where 
they were revealed [1]. The baseline distance A£ was 
counted from the place at CERN where the proton beam 
time-structure is being ascertained (it is the so-called 
Beam Current Transformer or BCT) to the origin of the 
OPERA detector at LNGS. The coordinates in the Uni- 
versal reference frame ETRF2000 [9] of the two locations 
were obtained by two steps: firstly the coordinates of 
ancillary benchmarks placed outside the two laborato- 
ries were determined by a GPS campaign and secondly 
the distance between these benchmarks and the BCT at 
CERN or the OPERA detector at LNGS was directly 
measured by geodetic survey. The value of A£ was then 
derived from the coordinates by usual Euclidean geome- 
try [10] yielding 731278.0 ± 0.2 m (deriving it from the 
spatial part of the metric as in (2) would add a negli- 
bible correction). Effects from geoid undulation, crust 
tides and continental drift (including earthquakes) do not 
significantly modify the above result. After an accurate 
synchronization between clocks of both laboratories [11], 
the baseline length was divided by the time of flight of 
neutrinos to obtain their average velocity. The result 
was larger than c by an amount [{v — c)/c] experiment = 
2.48 ± 0.28(stat) ± 0.30(syst) • 10"^. However, the mech- 
anism described in the paper is unable to explain this 
excess. This can be easily seen by adopting the simplify- 
ing hypothesis that neutrinos travel over the sphere of the 
Earth's surface and taking advantage of the fact that the 
angular sector of the Schwarzschild metric in standard 
coordinates is flat. Then A£ equals the Euclidean result 
r©A6' {A6 is the angular separation between CERN and 
LNGS [13]) while At is r©A6i/y^l - 2TO©/r©. Hence v/c 
is always 1. Eliminating the previous simplifying hypoth- 
esis (neutrinos actually traversed the Earth's crust fol- 
lowing the imaginary chord that joins CERN and LGNS) 
adds a negligible contribution of opposite sign (to verify 
this assertion an interior metric was derived at first or- 
der in TO0 assuming a planet Earth with uniform mass 



density). The inclusion of the spin of Earth (for ex- 
ample by using the Lense-Thirring metric) produces an 
even smaller contribution (while TO©/r© ~ 7 • 10""^*^, 
Je/r^ ~ 4 • 10~^^, J© being the Earth's angular mo- 
mentum). Only the rotation of the observer's laboratory 
at Gran Sasso induces a revealable effect but it affects the 
synchronization of clocks and therefore it has nothing to 
do with the topic described in this paper [12]. 

These considerations seem to imply that the obser- 
vation of a superluminal velocity for the neutrinos of 
OPERA must be likely ascribable to other reasons: ei- 
ther purely experimental oversights [14,15] or really new 
physics peeping out (see [16-36] for a partial list of the- 
oretical suggestions and criticisms). 

Consider now the lunar laser ranging experiment 
(LLR) [37,38] in which, among other ephemeris, the 
Earth-Moon distance is determined. It consists in send- 
ing a laser pulse from the Earth to the lunar surface 
where it is reflected (several manned missions in the past 
left corner reflectors on the; Moon's surface) and received 
back on Earth. Multiplying the time employed by the 
light in its round trip by the speed of light yields the 
Earth-Moon distance. However, assuming that the met- 
ric is well-approximated by the Schwarzschild one at lin- 
ear order in m^/r^, the mean speed of the light ray turns 
out to be larger than c and so, the distance should come 
out less than what it really is. To remedy this inconve- 
nient and obtain at least an order of magnitude of the 
necessary correction, we resort again to expressions (6). 

Establishing that the radial coordinate of the labora- 
tory is the Earth's radius r©, it has yet to find the radial 
coordinate on the Moon, n. This is achieved by invert- 
ing the formula for At in (6). After some algebra we get 
the equation 



exp 



At 
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(8) 



where ^ = ri/r©, p© = r©/(2TO©) and At is half of 
the time (measured by terrestrial clocks) spent during 
the round trip. Eq.(8) can be solved for ^ in terms 
of the Lambert T4^ -function [39]. However, on account 
of the fact that p© ^ 1, we can approximately set 
^ w 1 -|- AT/(r© — m®). Inserting it into the first of 
(6) leads to 
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(9) 



the correction being 2TO©Ar/r0 « 53 cm. It must be 
stressed that the above exercise has been presented in 
order to display an instance where the issue discussed in 
this paper yields a not insignificant contribution. With 
another name, this effect has surely been taken into ac- 
count by the LLR collaboration. Indeed, the time lapse 
during the round trip of the laser pulse was evaluated 
in barycentric reference frame coordinates [38,40] by in- 
cluding all general relativistic effects and leading to the 
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Shapiro formula [41]. Moreover, in the original calcula- 
tion of the LLR collaboration, the effects from all main 
Solar System bodies are included. 

In conclusion we have seen that mean velocities in gen- 
eral do not conform to well-known special relativity prin- 
ciples. In particular, albeit amazing, the average velocity 
at which a particle has travelled for a long time can be 
different from the instantaneous velocities that the par- 
ticle attained at every point along the trajectory, even 
when the latter was the same at all points. Indeed, spe- 
cial relativity adequately describes physics only locally 
while the average over large distances of a velocity must 
necessarily be introduced as a non local quantity. Uti- 
lizing these considerations for understanding the results 
from the OPERA collaboration, we conclude that the 
real anomaly is not that {v — c)/che positive but that it 
is a rather large number. 

The discussion about mean velocities can be straight- 
forwardly generalized to time dependent metrics as sen- 
sible definitions of proper time and length are also ad- 
missible on such metrics. An instance is the Friedmann- 
Robertson-Walker metric. 

For the Schwarzschild metric, the effect described in 
the paper is driven mainly by the curvature dependence 
of the proper time At. Indeed, note that as pi ap- 
proaches 1, the value of At diverges while that of A£ 
stays finite for any values of pi and p2- This property 
is not general. For exa mple, i n a conformally flat metric 
g^u = f{x)Vtiiy, At = ^/fix2)At while A£ = J dx^ f{x), 
which shows that any divergent behavior of the func- 
tion f{x) may influence both At and An exam- 
ple of a conformally flat metric is again the Friedmann- 
Robertson- Walker metric, although it is more usually 
presented in Gaussian coordinates. 

We have analyzed the case in which the observer stays 
at rest at one end of the particle's trajectory. But of 
course other experimental dispositions are possible; the 
observer staying in the middle of the particle's trajectory 
or even at a point outside it; the observer not at rest, etc. 
Also many possible deflnitions of mean velocity, other 
than the one used by OPERA (1), can be conceived. In 
all cases bizarre results should be carefully interpreted. 



It is a pleasure to thank Mihail Mintchev and Gian- 
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Note added: During the refereeing of the paper we be- 
came aware of the existence of Ref. [42] where a similar 
analysis is performed. However the study presented here 
is more general and detailed. 
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